Scaling law and stability for a noisy quantum system 
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We show that a scaling law exists for the near resonant dynamics of cold kicked atoms in the presence 
of a randomly fluctuating pulse amplitude. Analysis of a quasi-classical phase-space representation of 
the quantum system with noise allows a new scaling law to be deduced. The scaling law and associated 
stability are confirmed by comparison with quantum simulations and experimental data. 
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Coherent quantum phenomena may now be routinely 
observed in ultra-cold neutral atoms manipulated by light 
fields detuned from atomic resonance. The unprecedented 
control of atomic dynamics afforded by these atom-optical 
techniques has impacted a number of fields significantly in 
the last decade. In practical terms, the realization of cold- 
atom fountain atomic clocks and atom interferometers is 
very important for precision measurements and metrology 
in general [1]. Other promising applications include the 
manipulation of atoms in optical lattices I2D with possible 
applications to quantum computing 0]. 

Aside from such practical applications, atom optics has 
also offered the means to create ideal experimental imple- 
mentations of model systems, in particular, the quantum 
kicked rotor known in this realization as the atom optics 
kicked rotor. The system and its variants has been studied 
by a number of groups worldwide 0, EL Sill, HI] due to the 
ease of observing such quintessential quantum phenomena 
as dynamical localization [9] and dynamical quantum reso- 
nance [10]. Recent interest in the quantum resonance phe- 
nomenon comes not only from a fundamental perspective, 
but also from the useful features of the resonance behavior. 
For example it has been shown that the resonance peaks ex- 
hibit sub-Fourier resonance scaling opening the pos- 
sibility of faster than Fourier signal detection using the res- 
onance phenomenon [15]. Additionally, our work has great 
relevance to similar proposals for precision measurements 
of the atomic recoil frequency ifTTh . 

The cloud hanging over all planned implementations of 
quantum technologies, is that of decoherence Ill2ll - inter- 
action with environmental degrees of freedom which leads 
to irreversible loss of phase coherence in quantum systems. 
In atom-optics systems, decoherence typically arises due 
to spontaneous emission and timing and amplitude fluctu- 
ations in lasers. Typically, decoherence must be treated 
statistically, and its effect is only made plain by simulat- 
ing quantum master-equations. However, in the case of 
the quantum kicked rotor, some progress has been made in 
treating the response to spontaneous emission decoherence 
through a quasi-classical scaling theory lfl3ll . In this case 
the dynamics of kicked atoms near a fundamental quan- 
tum resonance, dependent ostensibly on four parameters 



(kick number, strength, period, and spontaneous emission 
rate) is reduced to a stationary function of two scaled time 
variables, with a closed analytical form. The presence of 
this scaling belies the fact that moderate noise typically de- 
stroys quantum correlations and it might be thought that 
the scaling function in the presence of spontaneous emis- 
sion is an isolated case where decoherence is analytically 
tractable. However, here we show that a scaling exists in 
the same system in the presence of amplitude fluctuations. 
Most remarkably, the fundamentally quantum decoherence 
process can be visualised with a classical phase space pic- 
ture here. The noise changes the topology of the phase 
space in a way that makes clear which parameter regimes 
will exhibit robustness to decoherence. 

It is important to note that amplitude noise induced de- 
struction of quantum correlations has been proven for non- 
quantum resonance conditions 111 411 . This naturally leads to 
the assumption that away from exact quantum resonance, 
amplitude fluctuations will rapidly induce quantum deco- 
herence. The contrary was proved by a recent experiment 
11511 . but the cause of the stability near quantum resonance 
has remained opaque. We derive in the following a thor- 
ough theoretical understanding of this robustness based on 
a semiclassical scaling approach. Our theory compares 
very well with measurements of near-resonant motion. 

Experimentally, we realize a kicked atom system with 
noise by overlapping an optical standing wave with a sam- 
ple of cold atoms and pulsing the potential periodically. 
The height of the potential can be controlled by adjust- 
ing the optical power transmitted through an acousto-optic 
modulator. The system with amplitude noise may be rep- 
resented by the Hamiltonian iflrjl 



H{?) = ^ + k cos(z) ^(1 + R s W/t - s) , (1) 



where p is the atomic momentum in units of 2hk L , z is the 
atomic position scaled by 2k L , t' is time, and t is the total 
number of kicks. Amplitude noise enters in the factors R s 
which are random numbers distributed uniformly on the 
interval [-L/2, +L/2], where L is a noise level between 
and 2. The scaled kicking period r is defined by the 
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FIG. 1: (a) Experimental measurements of quantum resonance 
peaks as a function of e for noise levels L — (*), L — 0.5 
(□), L = 1.0 (o), L = 1.5 (A) and L = 2.0 (o). (b) Points 
show experimentally measured peak energies (circles), whilst the 
dashed line shows the theoretical formula Eq. {4). Crosses show 
simulation energies, which exactly agree with the theory. Sample 
error bars are plotted from shot-to-shot measurements, not taking 
into account systematic uncertainties in the absolute value of k. 



equation r = 8u r T, where u r = hk 2 L /2M is the recoil 
energy. The kicking strength is proportional to the optical 
standing wave intensity, and its measured value was k ~ 
4.3 or k ~ 2.8 for the two separate sets of experimental 
data considered here. The kicking strength varied by about 
10% across the atomic sample. 

In our experiments, a sample of cold Cs atoms was pre- 
pared in a standard magneto optical trap (MOT) 111 511 . The 
atom ensemble had an initial width in momentum of up to 
a p /(2hk L ) ~ 8. They were released from the trap and ex- 
posed to either 5 or 20 periodic pulses of width 480 ns from 
an optical standing wave detuned by 0.5 GHz from atomic 
resonance. For the 20 kicks experiments (with k ~ 2.8) 
the presence of spontaneous emission at a rate of 2.5% 
per kick lead to a slight lifting and broadening of the reso- 
nance peaks. We corrected for the broadening by subtract- 
ing an additional small, empirically determined constant 
from the off-resonant energies in this case. Atoms were 
than allowed to evolve freely for 12 ms before applying 
the MOT beams and imaging the resultant fluorescence on 
a CCD camera. In this way, the momentum distribution of 
the atoms was calculated allowing comparison with theo- 
retical predictions. It has been shown that for pulse periods 
r equal to integer multiples of 2ir (so called fundamen- 
tal quantum resonances) a semi-classical map may be used 
to describe the quantum dynamics lfl3ll . We define a de- 
tuning e = t — 2tt£ which measures how far the pulse 
period is from the £th fundamental quantum resonance, 
and define new scaled momenta and position variables 
J s = \e\p s + ir£ + t(3, (where p s is the atomic momentum 
in units of 2hk L at kick s and (3 is the non-integer quasi- 
momentum) and 9 = z + 7r(l — sign(e))/2 mod(27r). 



Then the pseudo-classical standard map with amplitude 
fluctuations is (see Il3l[l5ll ) 

J s+ i = J s + |e|fc(l + R S ) sin(0 s+1 ) , 9 S+1 = 9 S + J s . 

(2) 

We now proceed to investigate how the mean en- 
ergy at exact quantum resonance is affected by am- 
plitude noise. To do this we need to find the aver- 
age over all amplitude noise realizations (and later ini- 
tial conditions 9 , J ) of the equation E t (k 1 , . . . ,k t ) = 



r(J t - Jo) 



E' s =o k *+i sin(0 o + s J Q ) 



where we have used an expansion given in fl!3ll . 
The noise average is given by (E t (ki, . . . ,k t )) R . = 

n*=i l S-l/2 dRj E t(h, h)- Since the series {k s = 
k(l + R s )} is a series of independent random variables, 
this expression simplifies greatly. Noting that (Rj) = 
and (RjRi) = 0, j / i, we need only retain the follow- 
ing terms of E t in the integrand: 



y^sin(6> + sJ ) 



s=0 



+ 



t-1 

s=0 



sin 



(3) 



We note in addition that (El=o ^ sin(0 o + s Jo)] 2 ) = 
12 S s =o sin 2 (6*o + sJ ), where we have used the fact that 



(l/L*) J^ L 2 /2 dRi . . . dR s ...dR t R 2 s = L 2 /12. Averag- 
ing over initial conditions (9 , Jo) gives, with 9 £ [0, 2ir] 
and J € [ir£, tt£ + r] corresponding to a uniform quasi- 
momentum distribution in the unit interval (see lfl3ll ): 
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(4) 



where we have used the fact that the averages over both 
terms in Eq. (([3]l) evaluate to t/2. (This result was also 
given in ref. Il7ll from a purely quantum argument). Fig. 
1 shows experimental data compared with simulation re- 
sults and Eq.[4] demonstrating good agreement between all 
three. Shot to shot errors were found not to vary with e and 
the given errorbars are estimates calculated from the stan- 
dard error over 10 energy measurements at a kicking period 
of 58/US. The discrepency between theory and experiment 
in the L = case is due to the difficulty in measuring the 
high momentum components, a problem which is amelio- 
rated by the addition of noise 



We now show how the scaling law introduced in II13H can 
be modified to take amplitude noise into account. We start 
with the pseudo-classical scaling function II 1 3fl 



(E t ,o) 



= R(t,k,e) « H(x) = 1 



$o(x) + — G(x), 

7TX 

(5) 

where x = tyjk |e| and (E t .o) = k 2 t/4: is the mean peak 
energy. The functions <T?o and G are evaluated numerically, 
and the reader is referred to ref. lfl3ll for details. 

For L > 0, we generally expect a loss of the scaling in 
all the variables e, k, t, L due to higher correlations in the 



FIG. 2: Phase space diagrams showing the effect of amplitude 
noise on the pseudoclassical map (2j. Left panel: without noise. 
Middle panel: noise level L = 1.5. Right panel: pendulum tra- 
jectories for various initial conditions. The separatrix is shown 
by a thick, dashed line. The grey shaded area shows a region of 
± ^AJ t 2 cs ^ for L = 1.5 about the separatrix, demonstrating the 
trajectories which lead to the correction in ©. 

evolution of the classical map neglected above when 
deriving (©. Remarkably, however, by observing the type 
of change in topology of the pseudo-classical phase space 
when increasing e, as depicted in Fig. [2] we can neverthe- 
less accurately estimate the change of energy growth in the 
presence of noise for small e (for which the semiclassical 
approach is valid for long experimental time evolutions). 
Noise is well known to enhance diffusion along nonlinear 
resonances in the first place [18;]. Therefore, we expect 
the major contribution of energy enhancement around the 
separatrix region of pseudo-classical phase space, which 
separates the two different topologies that give rise to the 
contributions G and <J> to the scaling function [13]. Since 
G describes bounded librating pendulum motion within the 
principal resonance zone, local changes of that motion due 
to noise will be small. The largest perturbation comes from 
classical trajectories moving close to the separatrix which 
is washed out due to the fluctuations of k (see Fig. |2). In 
this region, trajectories can actually perform rotating mo- 
tion now, where at L = they would still be bounded to 
the resonance. The increase of energy arising from those 
trajectories can be estimated by considering the area in 
phase space covered by them, as shown for L = 1.5 in 
Fig. 12 Since the width of the principal resonance is given 
by A J rcs ~ Ayke, the relative change in weight of rotat- 
ing orbits is given by 

1 (AJ r 2 cs ) - A J 2 S ^ L_ 
2vr AJ 2 S ~ Svr ' 

The noise- averaged standard deviation is (AJ^ S ) ~ (1 + 
L/4) A J r 2 cs by a simple integration. With this result we can 
now add the additional energy of rotating trajectories to the 
scaling function from (O, by adding a term L/(8tt)<&(x). 
Dividing now the true energies by the result at exact quan- 
tum resonance and L = 0, we finally arrive at the new 
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FIG. 3: The theoretical scaling function for L — (thick solid 
line), L = 0.5 (dotted line), L = 1.0 (dashed line), L = 1.5 
(dot-dashed) and L = 2.0 (thin line). Simulation data is also 
shown, rescaled by the factor \k 2 t, for L — 0.5 (□), L = 1.0 
(o), L — 1.5 (A) and L — 2.0 (o). Open symbols are produced 
for fixed k = 2.8, varying e G [10~ 3 ,0.1],t G [20, 150], while 
filled symbols represent data for randomly chosen values of k G 
[1, 10], e G [10 -3 , 0.1], t G [2, 150]. 

scaling function for finite noise 

= 1 + ^- - (1 - L/(8tc)) * (a:) + — G{x) . (7) 
12 irx 

Our derivation of (0 is thus analogous to the noise-free 
case, taking into account, however, the main contribution 
of heating due to noise. Higher-order correlations and heat- 
ing of the librating modes are neglected. We note that the 
principle changes to the phase space which give rise to this 
scaling are readily seen in Fig|2] In essentials, the scaling 
functiuon reduces a complicated quantum system which in- 
cludes decoherence to the dynamics of the pendulum. 

Inspection of Eq. (0 reveals some interesting features 
as seen in Figj3j Firstly, because <I>o saturates to 1 and 
G(x) is small for small values of x, the small x behavior is 
largely unchanged in the scaling function. Essentially, the 
zero-noise scaling function is merely displaced upwards for 
small t, k or e. Experimentally, this means that as long as 
x = t\Jk\ e| < 4 (e.g. take i = 20, A; = 0.1 and scan 
over e for any noise value), the resonance peak will not 
be broadened. This fact is important for proposed preci- 
sion experiments such as 111 ill where experimenters need 
to know how much tolerance the resonance width has to 
naturally occurring laser power fluctuations. Secondly, for 
large x the scaling function is significantly changed with 
the offset being much greater, corresponding to real broad- 
ening of the peak and reduction of peak visibility. 

A comparison of the theory with simulation results is 
shown in Fig. [3] It may be seen that the scaling function re- 



4 




1 10 1 

x 



FIG. 4: The theoretical scaling function of Eq. I0 (as shown 
in Fig. [3]) is compared with rescaled experimental data (as, e.g., 
from Fig.Q]i. Shown are data across more than one order of mag- 
nitude in the scaling variable x for (a) L — (*), L = 1.0 (<>), and 
L = 2.0 (o) and (b) L = 0.5 (□) and L = 1.5 (A). Theoretical 
curves are shown with the same line-styles as in Fig. [3] Note that 
for x < 5 the data comes from separate 5 kick experiments, and 
the error bars are the same size as the plotted points. Sample error 
bars are calculated as described in the caption of Fig.Q] Data for 
|e| > 0.15 is excluded since the pseudo-classical theory breaks 
down in this region. 



produces the broad shape of the quantum simulations over 
a large spectrum of parameters. Each point in Fig. [3]is ob- 
tained by averaging over 50,000 initial conditions, each of 
which is subject to kick-to-kick amplitude fluctuations. Al- 
though our statistics are good, there is still a non-negligible 
scatter in the simulation data which decreased systemat- 
ically when augmenting the number of initial conditions 
averaged to obtain the final energy. The experimental data 
from Fig.Q](a) and additional new data sets have been plot- 
ted in Fig. [4] The experimentally measured energies are 
obtained as an ensemble average over the total number of 
atoms and are rescaled by subtracting the mean initial en- 
ergy of the ensemble /4 and then dividing by the energy 
at the peak maximum for L = 0. The estimated error bars 
shown in Fig.[4]represent shot-to-shot fluctuations over dif- 
ferent noise realisations calculated as for Fig. [TJ 

In summary, we have derived and tested a generalized 
scaling function for the quantum resonance peaks in the 
presence of noise. The theory shows broad agreement with 
both quantum simulations and experimental results. Most 
importantly it illuminates new facts about the response of 
quantum resonance to noise - in particular, the stability of 
motion near to quantum resonance is revealed to be due 
to the unexpected persistence of scaling laws in the noisy 
system. Although the effect of amplitude noise is to mod- 
ify and even destroy quantum correlations, the effect near 
to quantum resonance can be understood precisely in terms 
of the noise-induced changes to the epsilon-classical phase 
space. Hence, quantum decoherence may be understood by 
a quasi-classical analysis in the system studied here. The 
robust nature of the scaling for small x allows us to predict 



parameter families of t, r and k for which noise will have a 
minimal effect on the quantum resonance, and surprisingly 
we find that for small enough x, the quantum resonance 
peak shape is entirely unaffected by noise (although a dis- 
placement in energy occurs). The exploration of quantum 
systems which exhibit resistance to noise is of great impor- 
tance for the future of quantum technologies. Analytical 
methods for determining the response of a system to noise 
and perturbations, as done here and in a different context in 
11911 . are valuable because they offer insights on stability of 
quantum motion which simulations cannot readily provide. 
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